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We calculate the full one-loop electroweak corrections to tri -boson production (ZZZ and WWZ) at 
the ILC. This is important to understand the Standard Model (SM) gauge quartic couplings which 
can be a window on the mechanism of spontaneous symmetry breaking. We find that even after 
subtracting the leading QED corrections, the electroweak corrections can still be large especially 
as the energy increases. 
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1. Introduction 

Due to its clean environment an e + e~ linear collider in the TeV range is an ideal machine 
to probe in detail and with precision the inner working of the electroweak structure, in particular 
the mechansim of symmetry breaking. From this perspective the study of e + e~ — > W + W~Z and 
e + e~ — > ZZZ may be very instructive and would play a role similar to e + e~ — > W + W~ at lower 
energies. Indeed it has been stressed that e + e~ — > W + W~Z and e + e~ — > ZZZ are prime processes 
for probing the quartic vector boson couplings [jl|]. In particular deviations from the gauge value 
in the quartic W + W ZZ and ZZZZ couplings that are accessible in these reactions might be the 
residual effect of physics intimately related to electroweak symmetry breaking. Since these effects 
can be small and subtle, knowing these cross sections with high precision is mandatory. This calls 
for theoretical predictions taking into account loop corrections. 

In this report we discuss the next-to-leading order (NLO) corrections to the processes e + e~ — > 
ZZZ and e + e~ — > W + W~Z at the future international linear collider (ILC). We also address some 
technical issues related to numerically stable evaluation of one-loop integrals, which is a challenge 
for one-loop multi-leg automatic calculations. 

2. Tree level 




Figure 1: Representative Born diagrams for e + e~ — > ZZZ and e + e~ — > W + W~Z. Diagrams (a) contribute 
to both processes while diagrams of type (b) contribute only to e + e~ — > W + W~Z. The first diagram of type 
(a) will be referred to as the Higgsstrahlung contribution. 

At leading order W + W~Z and ZZZ final states are produced through the diagrams shown in 
Fig. [I]. The important gauge couplings fermion-fermion-vector, yWW and ZWW, which also appear 
in the well-tested e + e~ — > VK + VK~process, have been measured and found in good agreement with 
the prediction of the SM [Q, B]. Both processes include the Higgsstrahlung contribution where the 
splitting H* — > VV occurs. This contribution is small and vanishes in the large Higgs mass limit. 
Since the precision electroweak data suggest a Higgs mass below the WW threshold, we restrict 
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our study to the region Mh < 160GeV. This means that the Higgsstrahlung contribution can not 
be resonant and therefore in our calculation no width is introduced. The important features at tree 
level are that the neutral quartic gauge coupling ZZZZ vanishes in the SM and the charged yZWW, 
ZZWW couplings occur in the e + e~ — > W + W~Z process. In practice, the coupling yZWW (and 
also yyWW) can be tested at lower energies via the process e + e~ — > W + W~y. 



3. NLO calculations 



Our calculations are done in the framework of the SM. The virtual corrections have been 
evaluated using a conventional Feynman-diagram based approach using standard techniques of 
tensor reduction. We use the packages FeynArts and FormCalc-6 . to generate all Feynman 
diagrams and helicity amplitude expressions [|J]. We also use SloopS to check the correctness 
of the amplitudes by checking non-linear gauge invariance (see ^ and references therein). The 
total number of diagrams in the 't Hooft-Feynman gauge is about 2700 including 109 pentagon 
diagrams for e + e~ — > W + W~Z and about 1800 including 64 pentagons for e + e~ — > ZZZ. This 
already shows that e + e~ —> W + W~Z with as many as 109 pentagons is more challenging than 
e + e~ — > ZZZ. Indeed getting stable results for all scalar and tensor (up to rank 4) box integrals in 
the process e + e~ — > W + W~Z is a highly nontrivial task. The five-point integrals are reduced to 
four-point integrals by using the method of Denner-Dittmaier |Q] which does not involve the Gram 
determinant, detG = det(2p ; - • pj) with pij being the external momenta, in the denominator. The 
four-point (and three-point) tensor integrals are in turn recursively reduced to scalar integrals by 
using Passarino-Veltman (PV) method. The problem with this method is that the numerical results 
become unstable when the Gram determinant is small. For instance, the result for a rank-4 box 
integral includes tensor coefficients of the form 

Diik, = WGf' {3A) 

where the numerator is a complicated function of internal masses and external momenta whose 
indices have been excluded for simplicity. In many cases the function N vanishes in the limit 
detG — > 0, leaving the tensor coefficients finite. In particular, this is usually true if the internal 
particles are massive. This non-trivial behavior of the numerator which is a linear combination 
of scalar integrals can be spoiled by numerical cancellation or inconsistent approximations (like 
small mass/momentum approximations) in calculation of the scalar integrals, leading to numerical 



instabilities in the right-hand side of Eq. ( |3.1| ) when detG becomes small. A good way to solve 
this problem is therefore using higher-precision arithmetic in the calculation of loop integrals when 
numerical cancellation occurs. 

It is important to notice that the PV method fails when detG is exactly zero. If this happens 
the Appoint function of rank M can be written as a combination of (N — 1) -point functions of rank 
M. This is called segmentation [{7|]. We have exploited this fact to avoid the small detG region by 
using segmentation if the following condition is met 

^% < 10-, (3.2, 



(2p 
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Figure 2: Dependence of O^Jj W Zy on the soft cutoff 8 S (Ey < 8 s i/s/2) in phase-space slicing. Only 
the non-singular part is shown, i.e. the IR singular ln(m^) terms are set to zero. The result using dipole 
subtraction is shown for comparison with the error given by the width of the band. 
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where p} nax is the maximum external mass of a box diagram. This extrapolation is used only for the 
four-point integrals and turns out to be good enough for the present calculations. We have compared 
this to the method of using higher precision arithmetic (quadruple precision for the Fortran 77 code) 
and obtained good agreement. 

In addition to the virtual corrections we also have to consider real photon emission, i.e. the 
processes e + e~ — > W + W~Zy and e + e~ — > ZZZj. The corresponding amplitudes are divergent in 
the soft and collinear limits. The soft singularities cancel against the ones in the virtual corrections 
while the collinear singularities are regularized by the physical electron mass. To extract the sin- 
gularities from the real corrections and combine them with the virtual contribution we apply both 
the dipole subtraction scheme and a phase space slicing method. The former is used to produce the 
final results since it yields smaller integration errors as shown in Fig. ||. Further details are given in 

0- 

It is well-known that the collinear QED correction related to initial state radiation in e + e~- 
processes is large. In order to see the effect of the weak corrections, one should separate this large 
QED correction from the full result. It means that we can define the weak correction as an infrared 
and collinear finite quantity. The definition we adopt in this paper is based on the dipole subtraction 
formalism. In this approach, the sum of the virtual and the so-called "endpoint" (see [^] for the 
definition) contributions satisfies the above conditions and can be chosen as a definition for the 
weak correction 



"weak — "virt "endpoint- (3-3) 

For the numerical results shown in the next section, we will make use of this definition. 

Before presenting our numerical results it is stressed that we have performed the calculation 
in at least two independent ways both for the virtual and the real corrections leading to two in- 
dependent numerical codes (one code is written in Fortran 77, the other in C++). A comparison 
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of both codes has shown full agreement at the level of the integrated cross sections as well as all 
the distributions that we have studied. Moreover, we have done detailed comparisons with other 



groups QIC , 11 1 and obtained good agreement. 



4. Numerical results 

To absorb large corrections to the electromagnetic coupling and universal corrections due to 
the isospin breaking effects we use 



a 



«G U — sin t>w , 



K 



«G M = ot(0)(l+Ar) 



(4.1) 



with denoting the Fermi constant and dw is the weak-mixing angle, at tree level. The explicit 
form of Ar at one-loop order together with all the input parameters are given in [§]. When calculat- 
ing the NLO corrections we have to subtract the one-loop Ar contribution to avoid double counting. 
Since the real photon corrections are proportional to 05 (0), we require the full NLO corrections to 
be of order 6{a? G a(0)). 
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Figure 3: The total cross section for e + e~ — > ZZZ ( top) and e + e~ — > W + W~Z (bottom) as a function of \fs 
for the Born, full ff(a) and genuine weak correction. The panels on the left show the Born, the full NLO 
and the weak correction. The panels on the right show the corresponding relative (to the Born) percentage 
corrections. 
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e + e 



ZZZ: 



As shown in Fig. ||] the tree-level cross section rises sharply once the threshold for production 
opens, reaches a peak of about 1.1 fb around a centre-of-mass energy of 600GeV before very 
slowly decreasing with a value of about 0.9 fb at 1 TeV. The full NLO corrections are quite large 
and negative around threshold, —35%, decreasing sharply to stabilise at a plateau around y/s = 
600 GeV with —16% correction. The sharp rise and negative corrections at low energies are easily 
understood. They are essentially due to initial state radiation (ISR) and the behaviour of the tree- 
level cross section. The photon radiation reduces the effective centre-of-mass energy and therefore 
explains what is observed in the figure. On the other hand the genuine weak corrections, in the 
scheme, are relatively small at threshold, —7%. They however increase steadily with a correction 
as large as —18% at ^/i = ITeV. 
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Figure 4: From top to bottom: distributions for the WW invariant mass and the rapidity of the WW system 
for e + e~ — > W + W~Z . The panels on the left show the Born, the full NLO and the weak correction. The 
panels on the right show the corresponding relative (to the Born) percentage corrections. 



e + e~ -> W+W-Z: 

Compared to ZZZ production, the cross section for e + e~ —> W + W Z is almost 2 orders of magni- 
tudes larger for the same centre-of-mass energy. For example at 500 GeV it is about 40 fb at tree 
level, compared to 1 fb for the e + e~ — > ZZZ cross section. For an anticipated luminosity of lab -1 , 
this means that the cross section should be known at the per-mil level. The behaviour of the total 
cross section as a function of energy resembles that of e + e~ — > ZZZ. It rises sharply once the 
threshold for production opens, reaches a peak before very slowly decreasing as shown in Fig. ||. 
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However as already discussed the value of the peak is much larger, ~ 50fb at NLO, moreover 
the peak is reached around y/s = 1 TeV, much higher than in ZZZ. This explains the bulk of the 
NLO corrections at lower energies which are dominated by the QED correction, large and negative 
around threshold and smaller at higher energies. As the energy increases the weak corrections get 
larger reaching about —18% at y/s = 1.5 TeV. This is similar to the result of ZZZ production and 
is consistent with the behavior of double-logarithmic Sudakov corrections. 

In Fig. |] we show the distributions in the WW invariant mass and the rapidity of the WW 
system. Due to photon radiation, in the full NLO corrections some large corrections do show up 
at the edges of phase space. However, even after subtraction of the QED corrections the weak 
corrections cannot be parameterized by an overall scale factor, for all the distributions that we have 
studied. 

5. Conclusions 

We have presented a calculation of the full next-to-leading order correction to the processes 
e + e~ — > W + W~Z and e + e~ — > ZZZ in the energy range of the international linear collider and 
for Higgs masses below the WW threshold. These processes would be the successor of e + e~ — > 
W + W~ in that they would measure the quartic couplings WWZZ and ZZZZ which could retain 
residual effects of the physics of electroweak symmetry breaking. With this in mind we have 
subtracted the QED corrections and studied the genuine weak corrections in the scheme. We 
find that the weak corrections can be large and increase with the energy. 
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